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Presence  of  shunt  current  in  fuel  cell  stacks  can  lead  to  corrosion  and  power  loss  problems.  The  objective 
of  this  paper  is  to  develop  a  stack-level  model  for  calculation  of  shunt  currents.  The  simulation  model 
was  based  on  electrical  circuit,  and  created  using  Simulink®  software.  The  Simulink®  results  were  val¬ 
idated  by  using  PSpice®  software  and  comparing  with  experimental  data  of  an  electrolyzer  stack.  The 
Simulink®  model  is  also  used  to  evaluate  the  effectiveness  of  a  protective  current  method  to  reduce  the 
shunt  current.  The  protective  current  method  was  found  to  effectively  reduce  shunt  current  of  a  100-cell 
stack.  The  Simulink®  model  also  shows  that  shunt  current  is  highest  at  end  cells  of  the  100-cell  stack 
examined,  suggesting  extra  care  be  applied  to  end  cells  for  corrosion  prevention.  Monte  Carlo  simulation 
was  performed  to  examine  sensitivity  of  variance  in  voltage,  manifold,  channel  and  cell  resistance  on 
calculated  shunt  current.  The  sensitivity  analysis  shows  that  calculated  shunt  current  is  most  sensitive  to 
variance  in  manifold  resistance  and  followed  by  variance  in  voltage,  electrolyte  and  channel  resistance. 
The  calculated  power  loss  due  to  shunt  current  is  within  1%  for  the  conditions  examined. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

A  fuel  cell  is  an  electrochemical  device  that  converts  chemical 
energy  stored  in  the  fuel  to  electrical  work.  Individual  fuel  cells  are 
often  connected  in  series  to  meet  voltage  requirements  of  intended 
applications.  When  liquid  electrolyte  is  used,  for  example,  in  phos¬ 
phoric  acid  fuel  cells  [1]  and  alkaline  fuel  cells  [2],  electrolyte  in 
individual  fuel  cells  is  shared  through  a  common  manifold,  and 
shunt  current  may  be  present  [1-9].  Shunt  current,  also  known  as 
bypass  current  [3],  current  leakage  [4],  or  parasitic  current  [6],  may 
be  present  through  the  common  electrolyte  due  to  presence  of  the 
conductive  paths  and  non-zero  electrical  field  potential  gradient. 
Similarly,  shunt  current  may  also  be  present  when  liquid  water  is 
used  as  coolant  in  polymer  electrolyte  fuel  cell  stacks.  Presence  of 
shunt  current  is  not  desired  because  it  can  lead  to  fuel  cell  corrosion 
and  power  loss  problems. 

Corrosion  might  occur  in  iron-water  systems  through  elec¬ 
trochemical  redox  reactions  [10].  Acid  environment  is  present  in 
phosphoric  acid  fuel  cells.  It  might  also  be  present  in  low  tem¬ 
perature  polymer  electrolyte  fuel  cells  when  flooding  occurs.  Base 
environment  is  present  in  alkaline  fuel  cells.  The  voltage-pH  equi¬ 
librium  diagrams  of  iron-water  systems  reported  in  [10]  suggest 
that  corrosion  may  occur  in  fuel  cell  stacks  when  stainless  steel  is 
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used  as  separator  or  manifold  materials.  Shunt  current  can  also  be 
present  when  impurities  or  ions  exist  in  cooling  water  that  circu¬ 
lates  through  bipolar  plates  of  the  fuel  cell  stack. 

Models  are  developed  to  predict  shunt  currents  [3-9],  for  exam¬ 
ple,  in  common  manifolds  and  in  bipolar  plates.  Shunt  current 
can  be  reduced  by  increasing  the  electrical  resistance  of  the  chan¬ 
nels  and  manifold.  This  can  be  accomplished  by  decreasing  the 
cross-sectional  area  and  increasing  the  length  in  the  channels  and 
manifold,  or  decreasing  the  resistance  in  the  cell.  The  cell  resis¬ 
tance  can  be  decreased  by  increasing  the  area  of  the  electrodes  and 
decreasing  the  distance  between  the  electrodes.  However,  they  will 
also  result  in  an  undesired  increase  of  the  pressure  drop  of  the  elec¬ 
trolyte  flow.  Another  method  to  reduce  the  shunt  current  is  to  apply 
protective  current  that  counters  the  electrical  voltage  gradient  in 
the  manifold  [11  ].  The  power  level  of  applied  protective  current  is 
usually  negligible  when  compared  to  fuel  cell  power  output. 

The  objective  of  this  paper  is  to  develop  a  stack-level  model  for 
calculation  of  the  shunt  current  distribution  in  fuel  cell  stacks.  An 
electrical  circuit  model  is  used  to  simulate  the  electrical  process  of 
liquid  electrolyte  fuel  cell  stacks.  The  calculation  is  compared  with 
experimental  data  of  an  electrolyzer  stack  reported  in  [11  ]. 

2.  Electrical  circuit 

Shunt  current  can  be  modeled  [  11  ]  using  electrical  circuit  shown 
in  Fig.  1 ,  where  Re  is  the  electrolyte  resistance  of  individual  cells,  Rc 
the  channel  resistance,  and  Rm  is  the  manifold  resistance.  A  small 
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amount  of  current,  or  shunt  current,  will  flow  through  the  channel 
and  manifold,  resulting  in  power  loss  and  possibly  corrosion  prob¬ 
lems  in  manifold.  To  apply  protective  current  to  eliminate  the  shunt 
current,  the  voltage  drop  through  the  manifold  (/<o^m)  must  be  the 
same  as  the  voltage  drop  through  the  electrolyte  (JRe).  When  this 
condition  is  met,  there  will  not  be  a  voltage  gradient  between  the 
cell  and  manifold.  Consequently,  no  shunt  current  will  flow  through 
the  channels.  Zahn  et  al.  [11]  derived  an  analytical  expression  for 
the  protective  current,  k0,  which  assumes  the  following  form: 


To  verify  Eq.  (1 ),  a  Simulink®  [12]  model  is  established  for  the  elec¬ 
trical  circuit  of  a  four-cell  stack  shown  in  Fig.  1.  The  Simulink® 
model  determines  the  current  flowing  through  the  cells,  channels, 
and  manifold  using  the  following  parameters:  Re,  Rc,  Rm,  V,  /,  and 
k0.  Simulation  is  run  using  a  Matlab®  function  [12],  and  analysis  is 
performed  to  obtain  the  following  quantities: 


n 


SCabs  =  JjCCQl 

1  =  1 

(2) 

PL  —  ^ideal  —  P\ 

(3) 

^ideal 

^ideal  = 

(4) 

n 

P  =  ^(EC),V-k0V 

(5) 

i=l 


Eq.  (2)  calculates  the  magnitude  of  the  sum  of  shunt  current  in  the 
channels,  in  which  i  denotes  the  channel  number  and  n  the  total 
number  of  the  channels.  When  protective  current  is  applied,  the 
calculated  value  of  SCabs  is  a  measure  of  the  effectiveness  of  Eq.  (1) 
to  eliminate  shunt  current.  Eq.  (3)  calculates  the  percent  power  loss 
due  to  shunt  current,  where  Picieai  is  the  ideal  power,  i.e.,  the  power 
level  when  there  is  no  shunt  current  present  that  is  calculated  by 
Eq.  (4).  Eq.  (5)  calculates  the  actual  power  delivered  from  the  stack, 


P,  i.e.,  the  cell  power  delivered  less  power  required  by  protective 
current. 

PSpice®  software  [13]  was  used  to  validate  Simulink®  results. 
The  model  parameters  used  in  the  calculation  were  provided  by 
R.J.  Bellows,  and  set  to  n  =  4,  Re  =  0.2  Rc  =  150  Rm  =  50  V=  1  V, 
and  /  =  1  A.  Electrical  current  calculated  by  Simulink®  and  PSpice® 
was  in  agreement  to  eighth  decimal  point  in  mA  without  protective 
current,  and  to  sixth  decimal  point  with  protective  current.  Detailed 
comparison  can  be  found  in  [14].  The  four-stack  model  is  modified 
for  calculation  of  the  10-cell  stack  electrolyzer  of  Zahn  et  al.  [11] 
and  expanded  for  calculation  of  a  100-cell  fuel  cell  stack. 

3.  Results  and  discussion 

The  shunt  current  results  of  the  Simulink®  stack  model  are 
compared  with  experimental  data  of  the  10-cell  electrolyzer  stack 
reported  in  [11  ].  The  electrolyzer  dimensions  were  provided  by  R.J. 
Bellows,  which  had  an  electrode  area  of  968  mm2  (1.5  in.2 ).  The  dis¬ 
tance  between  the  electrodes  is  estimated  to  be  31.75  mm  (1.25  in.). 
Based  on  electrolyte  conductivity  of  0.15  £2  cm-1  [14],  the  chan¬ 
nel,  manifold,  and  electrolyte  resistance  is  estimated  to  be  1000, 
150,  and  2.187  £2,  respectively.  The  cell  voltage  is  set  to  2.54  V  with 
applied  current  of  420  mA.  The  experimental  results  had  precision 
error  of  either  ±0.39  mA  or  ±1  mA,  depending  on  whether  the  range 
selector  of  the  current  meter  was  set  to  10  or  30  mA  full  scales. 
Experimental  data  with  no  protective  current  are  summarized  in 
Fig.  2,  in  which  error  bar  denotes  precision  error  of  ±1  mA.  Also 
shown  in  the  figure  are  the  Simulink®  results.  The  simulation  was  in 
general  agreement  with  experimental  data  and  bounded  by  exper¬ 
imental  uncertainty,  except  for  the  end  channels,  channels  1  and 
10.  The  calculated  shunt  current  in  channels  1  and  10  was  below 
experimental  data  by  31%.  The  simulation  also  shows  a  small  posi¬ 
tive  shunt  current  in  channel  5,  while  a  small  negative  shunt  current 
was  recorded  in  the  experiment.  The  recorded  negative  current  was 
within  experimental  uncertainty  of  ±1  mA. 

Simulation  results  of  experiments  with  protective  current  of 
17.5  mA  are  summarized  in  Fig.  3.  Simulation  shows  that  shunt  cur¬ 
rent  leaves  from  first  half  of  the  cells  (i.e.,  channels  1  through  5)  and 
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Experimental  Results  vs.  Simulation  Results  -  Shunt 
Current  w/o  Current  Protection 


Fig.  2.  Comparison  of  simulation  and  experimental  shunt  current  without  protec¬ 
tive  current:  Re  =  2.187  Q,  Rm  =  150  ft,  Rc  =  1000  £2,  V=  2.54V,  /  =  0.42  A;  Error  bars 
±1  mA. 

returns  to  second  half  of  the  stack  (cf.,  channels  6  through  10).  The 
magnitude  of  calculated  shunt  current  decreased  from  end  chan¬ 
nels  toward  central  channels  and  was  negative  symmetric  between 
channels  1  and  10,  2  and  9,  3  and  8,  4  and  7,  and  5  and  6.  Shunt 
current  was  very  difficult  to  measure  as  suggested  by  error  bars  in 
Fig.  3.  Experimental  data  lack  the  negative-symmetry  characteristic 
required  by  theoretical  considerations,  suggesting  further  investi¬ 
gations  are  needed  to  resolve  the  differences  between  simulation 
and  measurements. 

To  estimate  uncertainty  in  measurements  due  to  variance  in 
resistance  of  manifold,  channel  and  cell  electrolyte,  and  voltage, 
Monte  Carlo  simulation  was  performed.  Monte  Carlo  simulation 
with  sample  sizes  of 5000, 10,000, 15,000  and  20,000  were  reported 
in  Schaeffer  [14].  The  histogram  based  on  sample  size  of  15,000  is 
somewhat  different  from  that  based  on  sample  sizes  of  5000  or 
10,000,  but  it  is  nearly  the  same  as  that  based  on  sample  size  of 
20,000.  The  results  presented  herein  are  based  on  sample  size  of 
15,000,  unless  specified  otherwise.  In  simulation,  a  standard  devi¬ 
ation  of  5%  of  the  mean  value  was  applied  to  manifold  resistance, 
channel  resistance,  and  cell  voltage.  Setting  variance  to  5%  is  some¬ 
what  arbitrary;  however,  it  represents  the  95%  confidence  level 
that  was  used  to  determine  experimental  uncertainty.  Two  sets 
of  results  were  examined  in  Schaeffer  [14]:  one  was  for  the  con¬ 
ditions  without  protective  current  and  the  other  with.  The  results 
with  protective  current  will  be  emphasized  in  this  paper. 

The  histogram  of  the  results  with  variance  applied  to  manifold 
resistance  shows  that  channel  shunt  currents  varied  in  the  range 
of  ±0.3  to  ±0.4  mA.  There  also  existed  a  non-zero  finite  probability 
of  channel  current  to  reverse  its  direction  in  channels  5  and  6,  and 
to  a  lesser  degree,  in  channels  4  and  7.  When  variance  was  applied 


Fig.  4.  Relative  error  of  calculated  shunt  currents  as  compared  to  the  mean  values 
of  the  conditions  with  standard  deviation  of  5%  of  the  mean  applied  to  voltage  (dia¬ 
mond  symbol),  manifold  (triangle  symbol),  electrolyte  (square  symbol)  and  channel 
(circle  symbol)  resistance,  respectively:  the  shunt  current  in  mA  is  also  shown  (solid 
line);  protective  current  applied. 

to  channel  resistance,  the  channel  current  was  seen  to  vary  in  the 
range  of  ±0.03  mA  for  central  channels  and  ±0.1  or  ±0.15  mA  for 
end  channels  [14].  No  reversed  current  was  predicted.  When  vari¬ 
ance  was  applied  to  electrolyte,  the  resultant  variations  in  channel 
shunt  current  were  fairly  constant,  having  a  value  in  the  range  of 
±0.10  to  ±0.15  mA.  There  was  a  small  non-zero  finite  probability  of 
channel  current  flow  reversing  directions  in  channel  6.  Simulation 
with  protective  current  was  also  performed  to  evaluate  the  effects 
of  voltage  variations.  Simulation  shows  that  the  resultant  variations 
in  channel  current  are  fairly  constant,  at  a  level  of  approximately 
±0.3  mA.  A  non-zero  finite  probability  of  channel  current  revers¬ 
ing  the  direction  was  predicted  for  channels  5  and  6  and  channels 
4  and  7  (with  a  lower  probability  value).  Multiple  variances  were 
applied  in  Monte  Carlo  simulation  and  qualitatively  similar  results 
were  obtained  and  reported  in  [14]. 

Sensitivity  analysis  based  on  Monte  Carlo  simulation  was  per¬ 
formed  to  examine  how  variances  in  voltage  and  resistance  of 
manifold,  channel,  and  cell  electrolyte  impact  calculated  shunt 
current.  Sample  size  of  30,000  was  used.  Fig.  4  summarizes  the 
resultant  error  due  to  5%  variance  in  cell  voltage,  manifold  resis¬ 
tance,  electrolyte  resistance,  and  channel  resistance,  respectively. 
Variance  in  manifold  resistance  results  in  highest  shunt  current 
differences  across  all  channels,  followed  by  variance  in  voltage, 
electrolyte  resistance,  and  channel  resistance.  The  largest  differ¬ 
ences  occurred  in  channels  5  and  6  where  shunt  currents  were 
lowest  (as  shown  in  Fig.  4).  The  results  with  simultaneous  variations 
of  voltage,  manifold  resistance,  electrolyte  resistance,  and  channel 
resistance  are  summarized  in  Fig.  5.  Again,  manifold  variance  dom- 
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Fig.  5.  Relative  error  of  calculated  shunt  currents  as  compared  to  the  mean  values 
Fig.  3.  Comparison  of  simulation  and  experimental  shunt  current  with  protective  of  the  conditions  with  standard  deviation  of  5%  of  the  mean  applied  simultaneously 

current:  Re  =  2.187  £2,  Rm  =  150  £2,  Rc  =  1000^2,  V=  2.54  V,  7  =  0.42  A,  7P  =  0.0175  A.  to  voltage,  manifold,  electrolyte,  and  channel  resistance. 
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x  10  3  Chennel  Currents  without  Protective  Current 


Channel 


Fig.  6.  Shunt  current  of  a  100-cell  fuel  cell  stack. 

inates  the  shunt  current  calculation.  The  calculated  absolute  shunt 
current,  SCabs,  is  7.9  mA  and  corresponding  power  loss,  PL,  is  0.65%. 

To  test  the  effectiveness  of  the  protective  current  (Eq.  (1))  pro¬ 
posed  by  Zahn  et  al.  [11  ],  the  shunt  current  calculation  is  extended 
to  a  hypothetical  alkaline  fuel  cell  stack  of  100  cells.  The  resistance 
is  set  to  75  £2  for  each  manifold  segment,  500  £2  for  each  channel, 
and  0.25  £2  for  each  cell  [14],  which  assumes  separator  thickness 
of  1.2  cm,  manifold  length  of  9.0  cm,  channel  length  of  0.6  m,  and 
electrolyte  conductivity  of  194mScm-1.  To  calculate  the  wetted 
areas,  current  output  from  single  cells  is  set  to  1 A  and  current  den¬ 
sity  is  assumed  to  be  40mA cm-2.  The  resulting  wetted  areas  for 
separator,  channel  and  manifold  are  25, 0.61,  and  0.61  cm2,  respec¬ 
tively.  The  cell  voltage,  0.65  V,  was  estimated  based  on  EIS  model 
parameters  of  [15].  The  calculation  shows  that  the  overall  shunt 
current  is  23.9  mA  and  corresponding  power  lost  is  1.13%  of  deliv¬ 
ered  power.  The  shunt  current  distribution  across  the  100  cells  is 
summarized  in  Fig.  6,  showing  that  shunt  current  decreases  from 
end  cells  to  central  cells.  The  results  suggest  that  end  cells  of  the 
stack  would  experience  greater  risks  of  corrosion  problems.  When 
protective  current  of  12.0  mA  is  applied,  shunt  current  in  all  chan¬ 
nels  is  reduced  to  a  negligible  level,  10-15  A— lending  support  to 
robustness  of  Eq.  (1). 

4.  Summary  and  conclusions 

A  fuel  cell  stack  model  for  calculation  of  shunt  current  in  a  fuel 
cell  stack  was  created  using  Simulink®,  and  the  results  were  val¬ 
idated  by  PSpice®.  The  simulation  results  were  also  validated  by 
comparing  experimental  data  of  an  electrolyzer  stack,  which  was 
operated  with  or  without  the  presence  of  protective  current  [11]. 
When  no  protective  current  was  applied,  the  simulation  results 
compared  well  with  experimental  data.  The  agreement  is  gener¬ 
ally  within  experimental  uncertainty,  except  for  two  end  channels. 
When  protective  current  is  applied,  less  satisfactory  agreement  was 
obtained.  However,  discrepancies  were  generally  within  experi¬ 


mental  uncertainty,  except  for  three  end  channels.  The  protective 
current  proposed  by  Zahn  et  al.  [  11  ],  i.e.,  that  of  Eq.  ( 1 )  is  validated  by 
the  Simulink®-based  fuel  cell  stack  model  considered  in  this  paper. 
Monte  Carlo  simulation  was  performed  to  examine  how  variance  in 
voltage,  manifold,  channel,  and  cell  resistance  can  impact  the  stack 
performance.  Standard  deviation  at  5%  of  the  mean  was  applied. 
Major  findings  of  the  current  study  are 

(1)  The  protective  current  proposed  by  Zahn  et  al.  [11],  cf.,  Eq.  (1), 
can  effectively  reduce  shunt  current  in  the  100-cell  stack  under 
the  conditions  examined  in  this  paper. 

(2)  The  Simulink®  model  results  showed  that  highest  shunt  current 
was  present  at  end  cells  of  a  100-cell  stack,  suggesting  that  end 
cells  are  more  prone  to  corrosion  problems. 

(3)  The  Monte  Carlo  simulation-based  sensitivity  analysis  shows 
that  calculated  shunt  current  is  most  sensitive  to  variance  in 
manifold  resistance  and  followed  by  variance  in  voltage,  elec¬ 
trolyte,  and  channel  resistance. 

(4)  The  calculated  impact  of  shunt  current  on  power  loss  is  about 
1%  for  the  conditions  examined. 

(5)  Significant  differences  are  observed  between  measurements 
and  simulation  for  the  protective  current  case  shown  in  Fig.  3. 
Further  investigations  are  needed  to  resolve  the  differences. 
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